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, We study four-dimensional interacting fermions in a strong magnetic field, using the holographic 

Sakai-Sugimoto model of intersecting D4 and D8 branes in the deconfined, chiral-symmetric parallel 
phase. We find that as the magnetic field is varied, while staying in the parallel phase, the fermions 
exhibit a first-order phase transition in which their magnetization jumps discontinuously. Properties 
of this transition are consistent with a picture in which some of the fermions jump to the lowest 
. Landau level. Similarities to known magnetic phase transitions are discussed. 
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The study of QCD-like theories using a holographic approach has been a very fruitful line of investigation in the 
^ ■ last few years. Many qualitative features of QCD appear to be reproduced in these holographic models. In particular, 
I the Sakai-Sugimoto (SS) model [ll, based on the dynamics of D8-branes and D8-branes in the background of near- 
PsJ ■ horizon D4 branes, has yielded many interesting results. Most studies of the phase diagram Q, 13, d, 0, H S 

have focused on the confined phase or on the interplay between the chiral-broken phase and chiral-symmetric phase 
in the deconfined phase as a function of external parameters such as temperature, density, and electromagnetic fields. 

In this paper we will study relativistic interacting fermions at high density and strong magnetic field using this 
holographic description. This is done in the Sakai-Sugimoto model in the deconfined and chiral-symmetric phase. In 
' this phase the low-energy excitations are charged fermions. 
, As shown in ll|, [ij], when a magnetic field is turned on, the system behaves as if the charge is made out of two 
components, one which is represented holographically in the bulk by sources located at the horizon and the other by 
• '"j ■ smeared D4-branes inside the D8-brane along the radial direction. The first behaves under an applied electric field 
rS [ as charges in a dissipative neutral medium and the second behaves as if no dissipation is present. In this letter we 
■ show that, for a given temperature and density, there is a first-order phase transition as a function of the magnetic 
field for which there is a jump of the charges of the first kind to charges of the second kind along with a jump 
in the magnetization. A similar phenomenon in real materials where the magnetization grows rapidly is known as 
metamagnetism. The properties of the transition seem to indicate that it is consistent with a picture in which a 
phase transition of the fermions jumping to the lowest landau level has occurred. The fermions in the lowest landau 
level seem to form baryonic-like bound states which are singlets under the SU (N) and thus do not interact with the 
dissipative bath of gluons. 
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II. THE MODEL BACKGROUND 

We will consider the SS D4-D8 system in the parallel phase. The D4-D8 intersection has ND — 6 and so has as 
its lowest excitations only chiral fermions in four dimensions. In the phase where the D8-branes and Z?8-brane are 
parallel, chiral symmetry is unbroken and we have both left- and right-handed chiral fermions in the fundamental 
of SU{N). The fermions are charged under the U(l) symmetry on the D8-branes, and this interaction resembles 
electromagnetism but without dynamical photons. The system has two free parameters, L which is the distance 
between the D8-brane and the DS-brane, and the size i?4 of the circle in the direction. The system goes through 
various phases as a function of these parameters and external conditions such as temperature, density, and magnetic 
field. At low temperature the system is in the confined phase and at higher temperature T > ^® 
deconfined phase. In the deconfined phase at high enough temperature (which depends on L) the system will be in 
the chiral-symmetric phase. The phase transition we will find occurs at any temperature and L provided the density 
and magnetic fields are large enough. 

The metric and background fields of the near horizon limit of the D4-branes is given by 

ds^ = u^J{u){dx^f + d^^ + dxl)+u-"-(^J^^+u^dn'^ , 

e* = gsu'^'^ Fi^inia'f/^N.dVti (1) 

where 

/(-) = 1-4 (2) 

and u — Xfj, = and R — {ngsNcY^^y^. ut is related to the temperature, i.e. the inverse periodicity of the 
Euclidean time Xq , as ut = (47r/3)^i?^T^. In the holographic description the free energy is the action of the D8 
branes in the near-horizon limit of the D4-branes. 

Baryonic charge is identified with the charge under the U{1) living on the D8-branes. Adding charges to the system 
is accomplished by sourcing the component of the [/(I) gauge field. One can turn on a magnetic field in this U{1) 
as 2Tra F23 = h. At zero magnetic field the charges sit at the horizon, but at non-zero magnetic field, due to the 
Chern-Simons term in the world volume action of the D8-brane which holographically encodes the axial anomaly, Ai 
is turned on and induces a baryon charge proportional to h along the D8-brane [llj. In addition, an anomaly-driven 
axial current Ja = ^h^ is generated in the direction of the magnetic field. 

Generally, when a magnetic field is turned on, for a fixed L, there is a possibility that the system will favor a chiral- 
broken phase. However, we can change L at will and, for a given charge density (or chemical potential), delay breaking 
of chiral symmetry such that the magnetic phase transition occurs before the chiral phase transition. Whether this is 
true in QCD is hard to say. 

III. FOUR-DIMENSIONAL HOLOGRAPHIC FERMIONS 

The action of the D8 brane in this background includes both a Dirac-Born-Infeld (DEI) term and a Chern-Simons 
(CS) term which embodies the anomaly. Through the CS term the magnetic field couples the Ao and the Ai fields 
living on the D8 brane. We work with rescaled variables 



a 



M - —^A^ , h = 27ra'F23 (3) 
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•dbi = M duu''/^ ^(l - {a'^{u) f + f{u){a'M?) (l + 
Scs = -^-^ / {hao{u)a-^^{u) - haQ{u)ai{u)^ , (4) 



3 



where we have included both the 8-brane and anti-8-brane parts, with ao = oq and ai = — ai, so 



The integrated equations of motion are then given by 



= 6hai[u) + a 



v/i-KN)2 + /HK(u))2 



Vi-K(u))2 + /HKH)2 

where d is the total rescaled charge density.^ The boundary conditions are 00(1*7^) — and ai(oo) = 0, and there is 
no constant appearing in the equation of motion for oi due to regularity at the horizon. 

At non-zero magnetic field, some of the charge is found above the horizon and is represented by D4-branes smeared 
inside the D8-brane. Like ordinary baryons which holographically correspond to D4-branes localized on the D8- 
brane, these smeared D4-branes have Nc strings stretching between them and the D8-branes and so correspond 
holographically to bounds states of Nc quarks. The amount of charge located above a given u is given by — 3/iai(u), 
and the total charge in the form of smeared D4-branes is then just d» = —3hai{uT) ■ The rest, d^d^ — d+3hai{uT), 
remains at the horizon. 

After a little algebra, the equations of motion ^ can be put in the form 

'2 _ fjShai+d)^ 

"° ^ /u5 + fh^u^ + /(3/iai + d)2 - (3/iao)2 
r '2 ^ i^hapf 



A. Zero temperature 

At any non-zero temperature these equations can only be solved numerically. However, at zero temperature / = 1, 
so one can actually solve these equations analytically^ up to a transcendental equation which is solved numerically. Of 
course, at zero temperature the chiral-symmetric phase is never the preferred phase 13], but solving these equations 
in this case provides more control over the numerics and clarifies the results. As we will see, the numerical solution 
at low temperature is very close to the zero-temperature results. 

To solve the zero-tcmperature equations, one notices that dividing the two equations in ([7]) gives 

Oq 3hai + d , , 

a-i 6nao 

which can be integrated using the boundary conditions ai{oo) — and ao(oo) — fi. One then finds that 

(3/iai + df - (Shaof = d^ - • (9) 

By defining a new coordinate 

ih du , , 

, 10 



^ The baryon charge d is related to the physical density D by D = SZE^^rf. In these conventions, quarks (or, more generally, fermions) 

carry baryon charge 1 rather than 1/Nc- 
^ Note that at zero temperature one could in principle add a constant of integration to the once integrated a\ equation of motion which 

represents a current in the system, but this is not relevant to our case. 
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FIG. 1: A plot showing the three solutions to the integral equation H14|) for 2:00 with d — 1 and regulator e — 10"'^. The dashed 
red curve is Zoo while the solid blue curve is the regulated integral on the right hand side of (|14p . The values of Zoo at which the 
curves cross are solutions. The regulator is responsible for the plateau in the integral; as e is decreased, the plateau is pushed 
higher, and when e — > the largest solution is Zoo = 00. 

where Zoo = z{u = 00), the equations ([7]) reduce to 

d^ao ^ ai + — 
6h 

dzO-i = ao . (11) 
The solution to the equations of motion with the conditions 00(2 = 0) = and ai{z — Zqo) = is then 

d sinh z 



3/1 cosh Zr, 



From these equations we see that 



dcoshz d , 

«i = ^i. ^ ^ ■ (12) 

d/i cosh Zoo An 



n — — tanh Zoo 

«i(0) = ^(tZ^-^) ^ (13) 



3h \coshzo 



and thus, 







3h du , , 

(14) 



5 + h^u^ + 

cosn 2o. 



To find the actual sohition for z^{h, d), one needs to solve numerically. We regularize the integral^ by taking 
the lower limit of the integral defining z to be e and at the end sending e — > 0. For small d and h there are three 
solution for Zqo, but for large enough h only one remains which is Zoo ~ loge — > 00. Fig. [T] illustrates the case where 
there are three solutions. 

It should be noted that the appearance of three solution and the interesting phase structure which they imply is 
due to the use of the full DBI action ^ rather than just the Yang-Mills approximation. 

The free energy of the three solutions is shown in Fig. [2l At small h the solution Zoo ^ h is preferred, but at 
some finite h ~ 0.19 there is a jump to the branch parametrized by Zoo 00. This first-order phase transition is 
accompanied by jumps in both the magnetization and the chemical potential. 

As we can see, the preferred solution at large h has 



Having to insert by hand a regulator is a particular feature of T = 0. At non-zero temperature, uj- naturally cuts off integrals over u. 
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FIG. 2: Three solutions for (a) the free energy F and (b) the magnetization M at T = with d — 1. There is a first-order 
phase transition from the sohd blue z^a ~ h solution to the dashed red z^c oo solution at ft = 0.19. The dotted grey solution, 
which corresponds to the middle solution in Fig. [TJ connects the two stable solutions and is an unstable maximum of the free 
energy. 



We also find that ao ~ and o-i ~ ^jji £^re constant all the way up to w — > cx) where they jump to their asymptotic 
values ao = fi and oi — 0. Not only is all the charge smeared inside the D8-branes, = d, but it is in fact located 
right near the boundary at m = oo. 



B. Non-zero temperature 

We now turn to solving these equations in the more physical case of non-zero temperature. We will see that the 
properties seen for zero temperature persist. At high temperature there is only one solution, which means a single 
phase. As the temperature is lowered, there is a critical temperature Tc{d) below which one finds three solutions to 
the equations of motion for some range of h. One way to see this in the numerics is to plot ai{oo) as a function of 
ai(uT). At high temperature the graph (at all values of h) cross the horizontal axis only once. However, at lower 
temperature there is a range of h when there are three possible values of oi(ut) which gives ai(oo) = 0. This is shown 
in Fig. Ha). 
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FIG. 3: (a) Three solutions for Ai[oo) as a function of A]_{ut) for n = 1, T = 0.07, and h = 0.2. (b) The free energy F as a 
function of h for n — 1 and T = 0.09. There are two stable branches, the solid blue curve and the dashed red curve. The phase 
transition, where the blue and red curves cross, is at h — 0.235. 



Again, one of the solution is unstable and two are stable. The free energy of the two stable solutions as a function 
of h can be seen in Fig. [S^b). We see that there is a phase transition at some non-zero value of h. At a temperature 
slightly higher then the critical temperature both fi{h) and M{h) develop what looks like a wiggle; see Fig. ID As 
the temperature is lowered, we enter a regime where there is a phase transition; the chemical potential and the 
magnetization become discontinuous, as shown in Fig. ([5|). At zero temperature, the jump of the chemical potential 
is such that the axial current ja = ^hji attains its maximum value of d/2 and does not change any more as h is 
increased. This implies that in the large h phase the system is fully polarized. At non-zero temperature, ja jumps 
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close to the maximum value and asymptotes to it at large h. Similarly if we had turned on an axial chemical potential 
rather then a vector chemical potential, there would be a jump in the vector current in the direction of the magnetic 
field. 

M 




FIG. 4: (a) /i and (b) M are smooth functions of h for d = 1 and T — 0.102 which is just above the critical point at T = 0.101. 




FIG. 5: (a) /i and (b) M as functions of h for d = 1 and T — 0.09, below the critical point. There are now two branches of 
stable solutions, and the phase transition between them occurs at h — 0.235. 



To further describe the phase transition and to try to understand it better, we look at the value of the magnetic 
field at which the transition occurs as a function of the density and at fixed temperature. The results are plotted in 
Fig. [51^a) and closely fit 

h - _ (^IQ^ 

We also find that even at non-zero temperature the free energy of the phase at large h is well-approximated by 

As was expected, the behavior here is similar to the zero-temperature case but with finite temperature corrections. 
One qualitative difference is in the location of the charges in the high h phase. While almost all the charges are in 
the form of smeared D4-branes, they are no longer all near the boundary at u = oo but are instead spread across all 
values of u. 

In terms of the behavior as a function of the temperature, we find that for a given charge density d above a certain 
temperature there will be no phase transition for any h; thus, the line of first-order phase transitions ends at a critical 
point above which there is a smooth cross-over. The critical line in the T — d plane above which there is no phase 
transition is given in Fig. ^h) and behaves at large d as 

T-logd. (18) 

It appears that no matter how high the temperature is, one can find a large enough density for which a phase 
transition will occur at some strong magnetic field. Of course, at these parameters, for a given L, the preferred phase 
may be the chiral-broken phase. However, as explained before, we can always take L to be large enough that the 
chiral-symmetric phase is preferred. 
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C. Transport properties 

The conductivity of this system was computed in [l^] with the result 



^longitudinal — 



^ - ^rTT^^T 

where, as before, = —ihaiiuT) is the charge carried by the smeared D4 branes. These charges, located above the 
horizon in SU{Nc) singlet states, behave as if they do not interact with the dissipative bath. 

The phase transition discussed above is a transition in which changes discontinuously. Consequently, the trans- 
port properties also change discontinuously; the Hall conductivity grows and the longitudinal conductivity shrinks. 
Note that at zero temperature, in accordance with Lorcntz invariance, at any non-zero h the longitudinal conductiv- 
ity is zero and the Hall conductivity is ^d/h. Only at non-zero temperature do the different interactions with the 
dissipative bath result in a discontinuous conductivity. 

Because the D4-branes are essentially smeared baryons made of a delocalized collection of quarks, one might have 
thought that gluons would at some level be able to interact with the individual quarks, leading to dissipation. However, 
the thermal bath is holographically associated with the horizon, and since the smeared D4-branes are located above 
the horizon, their size in the boundary theory is much smaller than the wavelength of the thermal gluons. As a result, 
as seen by the gluons, the D4-branes are effectively localized SU{Nc) singlets with which they don't interact. 



IV. SOME INTERPRETATION 



To try to understand some of the physics of this transition, we explore the phase that is dominant at large h. At 
least some of its properties are similar to the properties of a free, massless fermion gas in four dimensions in its lowest 
Landau level. 

A massless fermion in four dimensions with charge e in a magnetic field B has the following energy spectrum: 

E = sJkl + 2eBn (20) 

where fci is the momentum in the direction of the magnetic field and n is a natural number labeling the Landau levels 
in the transverse directions. The degeneracy of each state is roughly B. If all fermions are in the lowest Landau level, 
the energy of the gas is given by 



(21) 
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and the density of fermions is 



giving a chemical potential and energy 



^/ ^ BkjYiax 7 (22) 



T,2 



.-i, E-^. (23) 

The lowest Landau level is not spin degenerate, thus when there's a jump to the lowest Landau level, there is an 
increase in the magnetization. We can estimate when the system crosses from the lowest to the first Landau level; this 
occurs when it is energetically favorable to increase n rather than add to the parallel momentum, meaning k'^^^^ ~ B, 
which corresponds to 

B ~ n)'^ . (24) 

Many of the properties of the large h phase agree with those expected from free fermions in the lowest Landau level 
in four dimensions. The behavior of the free energy and chemical potential (|15p match the free fermion result (I23p . 
and phase transition (|16p resembles that of a transition into the lowest Landau level 

On the supergravity side, in the large h phase almost all the baryonic charge^ is represented by smeared D4-branes 
that crept up the D8-brane. Recently, there have been other indications that smeared D4-branes seem to be related 



to fermions in the lowest Landau level. For example, in 1^ it was argued that in the confined, large ft., pion-gradient 



phase the smeared D4-branes also behave as if they are free fermions and reproduced the thermodynamics of the 
lowest Landau level. 



As shown in n| , in the chiral-symmetric phase a baryon chemical potential in a background magnetic field generates 



an axial current. This ano maly -driven current is also present in a gas of free fermions and is due only to the fermions 



in the lowest Landau level 15| . In the supergravity description, since there is no constant of integration in equation 
([6]), the fivc-dimcnsional current is located above the horizon; consequently, the axial current is due to the smeared 
D4-branes rather than the charges at the horizon. 

Furthermore, the smeared D4-branes holographically correspond to bound states of Nc fermions, making up what 
look like smeared baryons. This phenomena may be related to the fact that in the lowest Landau level fermions are 
confined to move only in one spatial dimension, and thus any attractive potential will make them bind, even one 
which would not bind them in three spatial dimensions. Such an occurrence has been demonstrated to happen for a 
different situation related to the magnetically-induced chiral-symmetry breaking [l6| . 

Because the Sakai-Sugimoto model provides a robust holographic model of QCD, it would be interesting to see if 
the phenomena described here could be found, for example, at RHIC in a strong magnetic field. We see the onset 
of anomaly-driven behavior as soon as a magnetic field is turned on, but this may be an artifact of exactly massless 
quarks. In real QCD, a very strong magnetic field may be required before any of these phenomena become manifest. 

There are a few non-relativistic fermionic systems which resemble in some way the relativistic system discussed here. 
In particular, several condensed matter systems also exhibit first-order quantum phase transitions with discontinuous 



magnetizations. In the case of metamagnetism 17[, the magnetization, which comes from spin ordering, jumps due to 
a large increase in the density of states near the Fermi surface. This effect sometimes occurs in paramagnetic materials 
with itinerant electrons in their outer shells. The magnetization, as a function of the external magnetic field, can 
have a large rise and sometimes a discontinuity due to a rapid transition of many fermion spins from anti-parallel to 
the magnetic field to parallel. 

Another phenomenon that can cause a jump in the magnetization is based on the de Haas- van Alphen (dHvA) effect 
which is due to orbital angular momentum. At low temperature free fermions in four dimensions develop oscillations 
in their magnetization and chemical potential, among other quantities, as a function of the inverse magnetic field. ^ 
These oscillations occur when the the Fermi surface sweeps through the Landau levels. 

Because the charges feel both the background magnetic field and field produced by the other charges, the magneti- 
zation M of the material is a function of the magnetic inductance _B, which itself is a function of both H and M . The 



* At T = this will be all the baryonic charge. 

^ For a recent computation of these oscillations for massless fermions see [T^ . 
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function M{B) = M{H, M) may therefore have more then one solution for M for a given H if the dHvA oscillations 
become large |19l . [2o| . In this case, a first-order phase transition may occur when the system transitions from one 
Landau level to the next. This phase transition, which happens in real metals only at very low temperatures, has 
been observed mainly in diamagnetic metals and is therefore called a diamagnetic phase transition.^ Depending on 
the geometry of the sample, the phase transition may lead either to a homogenous phase or to an inhomogeneous 
phase, with multiple Condon domains of different magnetisms. 

The phase transition wc find in the SS model appears to have some resemblance to metamagnetism but also some 
differences. We start with a paramagnetic state and see a jump in the magnetization at a particular value of the 
magnetic field. However, after the jump, the magnetization decreases, while in materials exhibiting metamagnetism, 
the phase after the transition is ferromagnetic. In addition, while we find evidence the large magnetic field phase in 
the SS model is associated with the lowest Landau level, metamagnetism, which occurs at more modest magnetic 
fields, is not related to lowest Landau level physics. 

The diamagnetic phase transition is associated with transitions between Landau levels, and it occurs each time the 
Fermi surface crosses a Landau level; we, however, see only one phase transition rather than a whole series. Indeed, 
in our system we do not observe the dHvA oscillations at all, even at zero temperature, which may indicate that 
there is no well-defined Fermi surface.^ Furthermore, in the SS model, we have only a non-dynamical background 
electromagnetic field, thus B = H . Hence, this is c^uite different than the diamagnetic phase transition. 

We are left then without a truly compelling interpretation of our results. In addition to serving as a holographic 
model of QCD, the SS model describes a system of strongly-coupled relativistic fermions, which to some degree 
reproduces phenomena seen in some condensed matter systems. There seems to be some analogy with metamagnetism, 
but this analogy is incomplete. 
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